Astronomy & Astrophysics manuscript no. model7s 


February 2, 2008 


(DOI: will be inserted by hand later) 





Toward Understanding Environmental Effects in SDSS Clusters 



o 
o 

(N 

> 
O 



> 

in 
m 
(N 



O 



o 



X 



J. Einasto\ E. Tago\ M. Einasto\ E. Saar', I. Suhhonenko^, P. Heinamaki^'^ G. Hutsi'-^^& D. L. Tucker'* 

' Tartu Observatory, EE-61602 Toravere, Estonia 

^ Tuorla Observatory, Vaisalantie 20, Piikkio, Finland 

^ Max-Planck Institut fiir Astrophysik, Karl-Schwarzschild-Str 1, 86740 Garching, ERG 
Eermi National Accelerator Laboratory, MS 127, PO Box 500, Batavia, IL 60510, USA 

Received 09.1 1.2004, Accepted 

Abstract. We find clusters and superclusters of galaxies using the Data Release 1 of the Sloan Digital Sky Survey. We calculate 
a low-resolution density field with a smoothing length of 10 Mpc to characterise the density of the cluster environment. 
We determine the luminosity function of clusters, and investigate properties of clusters in various environments. We find that 
clusters in a high-density environment have a luminosity a factor of ~ 5 higher than in a low-density environment. We also study 
clusters and superclusters in numerical simulations. Simulated clusters in a high-density environment are also more massive 
than those in a low-density environment. Comparison of the density distribution at various epochs in simulations shows that in 
large low-density regions (voids) dynamical evolution is very slow and stops at an early epoch. In contrast, in large regions of 
higher density (superclusters) dynamical evolution starts early and continues until the present; here particles cluster early, and 
by merging of smaller groups very rich systems of galaxies form. 

Key words, cosmology: large-scale structure of the Universe - clusters of galaxies; cosmology: large-scale structure of the 
Universe - Galaxies; clusters: general; cosmology: simulations; cosmology: evolution 



1. Introduction 

Clusters and superclusters of galaxies are the basic building 
blocks of the Universe on cosmological scales. The first cat- 
alogues of clusters of galaxies by Abell (IT9581 [T989t and 
Zwicky et al. (1961-681 were constructed by visual inspection 
of the Palomar Observatory Sky Survey plates. Modern sur- 
veys of galaxies, such as the Las Campanas Redshift Survey 
(LCRS), the Sloan Digital Sky Survey (SDSS) and the two- 
degree-field (2dF) Galaxy Redshift Survey, enable us to define 
groups, clusters and superclusters of galaxies and to investigate 
properties of these galaxy systems in various large-scale envi- 
ronments. 

Studies of the dependence of properties of galaxy sys- 
tems on the density of the large-scale environment have been 
made by Einasto et al. 2003a 2003b^ 2pp3c 2003d (hereafter 
E03a, E03b, E03c and E03d, respectively), using the Early 
Data Release (Stoughton et al. l2002> of the SDSS and the Las 
Campanas Redshift Survey. These studies demonstrated the 
presence of environmental efifects in clusters - clusters in a 
high-density environment are richer and larger than in a low- 
density environment. 

The present study has three goals. First of all, we shall 
check the results obtained by E03a and E03b using a more 
accurate definition of groups and clusters. E03a and E03b 
found clusters of galaxies as density enhancements in the high- 
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resolution 2-dimensional density field. Such a definition has its 
restrictions, as in some cases clusters may overlap in the pro- 
jection direction, and cluster properties may be distorted. In 
contrast to previous studies we shall now define groups and 
clusters in the conventional way using the full 3-dimensional 
data on the distribution of galaxies. In this analysis we shall 
use galaxy samples of the Data Release 1 of the Sloan Digital 
Sky Survey (DRl of SDSS, Abazajian et al. 2003), and shall 
investigate properties of these systems in relation to the large- 
scale environment, from rich superclusters to poor filaments of 
loose groups in voids. 

The second goal of our study is to compare properties of 
clusters and superclusters with properties of similar systems 
found in N-body simulations of structure evolution. In partic- 
ular, we are interested in the relationship of cluster proper- 
ties and their large-scale environment. We define DM-haloes 
in simulations in the same way as groups and clusters were 
defined in real galaxy samples, and superclusters as large over- 
density regions of the smoothed density field. This comparison 
of real and simulated cluster properties complements a similar 
study by Einasto et al. (2004b I where a different method was 
used to characterise the density of the environment. 

The third and ultimate goal of the present study is to try to 
find an explanation for the environmental dependence of clus- 
ter luminosities. For this purpose we shall compare the evolu- 
tion of groups and clusters in high- and low-density regions. 
Also we shall compare the distribution of particles located 
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in systems of various richness in high- and low-density envi- 
ronments. This comparison shall be done for various epochs, 
which gives us the possibility to follow the evolution in regions 
of various global density. 

In the next Section we describe the SDSS DRl sample of 
galaxies and the method used to find groups/clusters of galax- 
ies. Here we describe also the N-body simulations of the struc- 
ture evolution used to compare observations with models. In 
Section 3 we describe the density field of the SDSS DRl, and 
properties of clusters of the SDSS DRl. In Section 4 we com- 
pare the properties of observed clusters with those of similar 
objects found in simulations. In Section 5 we follow the evo- 
lution of high- and low-density regions in an attempt to under- 
stand the mechanism behind the environmental effects in clus- 
ter and galaxy luminosities. We discuss our results and com- 
pare them with previous studies in Section 6. The last Section 
brings our conclusions. High-resolution colored figures of the 
SDSS DRl density fields are available at the web-site of Tartu 
Observatory (http://www.aai.ee/~einastoi. Preliminary 
results of this study were reported at the conference on the Zone 
of Avoidance by Einasto et al. ( 12004at . 

2. Data 

2.1. SDSS DR1 data 

The SDSS Data Release 1 consists of two slices of about 2.5 
degrees thick and 65-105 degrees wide, centered on the celes- 
tial equator, and of several regions at higher declinations. In the 
present study we have used only the equatorial slices. We ex- 
tracted the Northern and Southern slice samples from the full 
DRl sample, using the following criteria: the redshift interval 
1000 < cz < 60000 km s"' , the Petrosian r*-magnitude interval 
13.0 < r* < 17.7, the right ascension and declination intervals 
145° < RA < 250.0° and -1.25° < DEC < 1.25° for the 
Northern shce, and 350° < RA < 55.0° and -1.25° < DEC < 
1.25° for the Southern slice. The number of galaxies extracted 
(Ng-ii) and the length of the slice in the right ascension (ARA) 
are given in Tabled 

The SDSS data reduction procedure consists of several 
steps: (1) calculation of the distance, the absolute magnitude, 
and the weight factor for each galaxy of the sample; (2) finding 
groups/clusters of galaxies using the friends-of-friends algo- 
rithm; (3) calculation of the density field using an appropriate 
kernel and a chosen smoothing length. When calculating lumi- 
nosities of galaxies we regard every galaxy as a visible member 
of a density enhancement (group or cluster) within the visible 
range of absolute magnitudes. Mi and M2, corresponding to the 
observational window of apparent magnitudes at the distance of 
the galaxy. This assumption is based on observations of nearby 
galaxies, which indicate that practically all field galaxies be- 
long to poor groups like our own Galaxy, where one bright 
galaxy is surrounded by a number of faint satellites. Using this 
assumption, we find groups and clusters with their haloes, as 
single giant galaxies with their companions, or groups/clusters 
with their faint members. Further, we assume that the luminos- 
ity function derived for a representative volume can be applied 
also for individual groups and galaxies. 



The calculation of the distances, absolute magnitudes and 
weight factors of galaxies was described in detail in E03a. 
When calculating total luminosities of galaxies on the basis of 
their observed luminosities we used the Schechter ( 197^ func- 
tion with three sets of parameters. One set is based on the SDSS 
luminosity function by Blanton et al. (2001 1 and is denoted B, 
the other two sets on the SDSS luminosity function found in 
E03a, denoted El and E2. The respective values of the charac- 
teristic luminosity M* and the shape parameter a are given in 
Tabled 

In Fig. [Owe show the luminous-density weights as a func- 
tion of distance. The weights by Blanton are rather large at 
large distances; the weights by Einasto El and E2 are lower. 
The difference between the sets El and E2 is due to the fact 
that selection effects influence the estimated total luminosities 
of clusters and superclusters in different ways. The weights E2 
have been derived with the aim to get the correct total lumi- 
nosity of the sample as a whole at a given distance from the 
observer; this set yields luminosities of superclusters indepen- 
dent of distance. However, in this case the visible clusters have 
to include also luminosities of the invisible clusters, and lumi- 
nosities of individual clusters are too high at large distances. 
To avoid this distortion of cluster luminosities we have used 
the weights of the set El; this set yields for clusters the esti- 
mated total luminosities, which are statistically independent of 
their distance (for details see E03a). The right panel of Fig. [2 
shows the observed and total luminosities of galaxies at various 
distances. 

The next step is the search for groups and clusters of galax- 
ies. Here we used the conventional friends-of-friends algorithm 
by Zeldovich, Einasto & Shandarin (19821 hereafter ZES). 
Another algorithm was suggested by Huchra & Geller ( 119821 
hereafter HG). These algorithms are essentially identical with 
one difference: ZES used a constant search radius to find neigh- 
bours whereas HG applied a variable search radius, depend- 
ing on the volume density of galaxies at a particular distance 
from the observer. We compared SDSS cluster catalogues ob- 
tained by both algorithms, and found that mean virial radii of 
groups/clusters are practically constant for the constant search 
radius, and increase with distance for the variable search ra- 
dius (for a comparison of group radii for both algorithms see 
Einasto et al. 2004a 1. In the following analysis we have used 
only the group/cluster catalogue found with a constant search 
radius. The number of groups/clusters found for both equatorial 
slices is given in TableQ] 

2.2. N-body models 

We used a flat cosmological model with the parameters de- 
rived from a joint analysis of the WMAP microwave back- 
ground experiment and SDSS data by Tegmark et al. (I2004> 
(see also Bennett et al. 2003 1. We calculated three models 
with cube sizes L = 100 and 200 h^^ Mpc. The smaller 
cube was calculated for a 128^ mesh and the same number of 
DM particles, and the two larger cubes for a 256^ mesh and 
particles; we designate these models as MlOO, M200A, and 
M200B, respectively. The cosmological parameters of mod- 
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Fig. 1. The left panel shows the weights of galaxies used to correct for invisible galaxies outside the observational window. Black 
symbols show the weights for the Einasto set El of the Schechter parameters (used to find total luminosities of clusters); light 
gray and dark gray symbols show the weights for the Blanton B and the Einasto set E2 of the Schechter function parameters, 
respectively (the weights calculated with the set E2 were used to find the low-resolution density field). In the right panel we plot 
luminosities of galaxies: dark gray symbols show the observed luminosities, light gray and black symbols show the expected 
total luminosities, obtained using the Blanton B and the Einasto El sets of the Schechter function parameters. 

Table 2. Data on N-body DM-haloes and superclusters 



Sample L[Mpc/h] Q„, as h N„ M„lW°h-' M^] N„„„ Nu,,o 

MlOO 100 0.30 0.70 0.80 0.70 128^ 3.98 8 2459 

M200A 200 0.30 0.70 0.80 0.70 256^ 3.98 

M200B 200 0.27 0.73 0.84 0.71 256^* 3.59 20 12306 



els are given in Table |2l here Q,„ is the matter density (dark 
plus baryonic matter), Oa is the dark energy density (all in 
units of the critical cosmological density), erg is the present 
density fluctuation parameter, and Mp is the mass of a sin- 
gle particle. Here and elsewhere h is the present-day dimen- 
sionless Hubble constant in units of 100 km s ' Mpc For 
the models MlOO and M200A the initial power spectrum 
was taken using the approximation formula given by Klypin 
et al. I I1993I I. For the model M200B the initial power spec- 
trum was generated using the COSMICS code by Bertschinger 
l |http : //arcturus -itiit . edu/cosmics); here we accepted 
the baryonic matter density Q.h = 0.044. 

In simulations we used the Multi Level Adaptive Particle 
Mesh (MLAPM) code by Knebe et al. <2001l l. This code uses 
an adaptive mesh technique in regions where the density ex- 
ceeds a fixed threshold. The DM-haloes were found using the 
conventional FoF algorithm with a constant search radius for 



haloes of density contrasts 6n/n — 80, 125 and 411; these cor- 
respond to neighbourhood radii b - 0.23, 0.20 and 0.134 in 
units of the mean particle separation, respectively. The density 
contrast 6n/n = 80 coincides with that used by Tucker et al. 
(2000) in the search of loose groups, 6n/n = 411 describes 
virialized haloes in our accepted "concordance" model (see 
Peacock 1999 1, and the intermediate value of the neighbour- 
hood radius 0.2 was advocated by Jenkins et al. <200H . The 
DM-haloes in model M200B, used subsequently in our analy- 
sis, were found using the neighbourhood radius b - 0.23. In 
Table 12] N,„in is the minimum number of particles in the DM- 
haloes, and Nhaio is the number of haloes found. 

In small DM-haloes some particles have rather large veloc- 
ities relative to the rest of the halo; these particles evidently 
do not belong to the virialized part of the halo. To avoid the 
inclusion of unbound objects we should apply the virial the- 
ory condition Er - Ekm/\Epot\ < 0.5 (here Epot is the potential 
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energy and Ej^in the kinetic energy of a group). However, in 
groups with too high kinetic energy only a small fraction of 
particles are responsible for this effect; thus by excluding all 
these groups we would reduce the number of groups too much. 
To reduce statistically this effect we applied in model M200B a 
more modest criterion, E,- < 0.8. The model M200A was used 
only to investigate the evolution of populations of particles of 
various local and global density, so for this model individual 
DM-haloes were not found. 

3. SDSS DR1 clusters 

3.1. The density field of the SDSS DR1 

The SDSS DRl equatorial slices are very thin, thus 
3-dimensional and 2-dimensional density fields are very similar 
to each other Taking this into account we calculated only the 
2-dimensional luminosity density fields for observational sam- 
ples. As in E03a and E03b we calculated the high-resolution 
luminosity density field using Gaussian smoothing with a rms 
scale of 0.8 /; ' Mpc, and the low-resolution field with a rms 
scale of 10 /; ' Mpc. The high-resolution field was found using 
the Schechter parameters of the set El, the low-resolution field 
with the parameter set E2. 

The low-resolution field yields information on large over- 
density regions. This field was used to define superclusters as 
connected over-density regions. As in E03a we used density 
thresholds of 1.8-2.1 to find superclusters. At lower thresholds 
superclusters start to merge into percolating systems, and thus 
the definition of superclusters as the largest but still isolated 
high-density regions is violated. For higher thresholds the num- 
ber of superclusters decreases rapidly, since many of them have 
lower peak density. The catalogue of superclusters found using 
the SDSS DRl with the Schechter parameters of the set E2 is 
rather close to the catalogue in E03a that used the SDSS EDR, 
except that the total luminosities of superclusters vary a bit due 
to the use of more complete data. 

3.2. Properties of the SDSS clusters in various 
environments 

Fig-Elshows the luminosities of groups/clusters at different dis- 
tances from the observer We see that there exists a well-defined 
lower limit of cluster luminosities at larger distances. The limit 
is linear in the log L - d plot. Such behavior is expected as at 
large distances an increasing fraction of clusters do not con- 
tain any galaxies bright enough to fall into the observational 
window of absolute magnitudes. Mi . . . M2. The limit is lower 
for groups containing only one galaxy in the visibiUty window; 
these groups are systems like our Local group with one bright 
galaxy surrounded by faint companions. The difference in the 
low-luminosity limits of groups containing at least one or two 
galaxies in the visibility window is by a factor of two as ex- 
pected (this factor corresponds to the case when two galaxies 
in the visibility window have the same luminosity). 

Let us compare now properties of groups/clusters in var- 
ious environments. We shall use the density found with the 
10 /i ' Mpc smoothing as the global density in the superclus- 
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Fig. 2. The luminosities of SDSS Northern groups/clusters at 
different distances, corrected for galaxies outside the visibility 
window. Grey symbols stand for groups with at least two vis- 
ible galaxies, black symbols - for groups containing only one 
galaxy in the visibility window. The luminosity dependence of 
galaxies of the SDSS Southern slice is very close to that of the 
Northern slice. 
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Fig. 3. The luminosities of the SDSS DRl Northern slice clus- 
ters as a function of the environmental density, found by 
Gaussian smoothing of the luminous density field with the rms 
scale 10 /z ' Mpc. Data for the Southern slice are very similar. 

ter environment of clusters. The environmental effect is shown 
in Fig.|3l There is a correlation between the luminosity of the 
most luminous clusters and their environmental density. The 
most luminous clusters in high-density regions have a luminos- 
ity a factor of about 5 higher than the most luminous clusters in 
low-density regions. Using a different definition of the large- 
scale environment, a similar effect was found in the vicinity 
rich clusters of galaxies by E03c and E03d, and in the vicinity 
of massive DM-halos by E04b. 

3.3. Cluster luminosity functions 

Fig. 0] shows the integrated luminosity functions of the 
groups/clusters of the SDSS DRl Northern and Southern sam- 
ples. The absence of low-luminosity clusters at large distances 
has been taken into account by the standard V^^. weighting 
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Fig. 4. Left: the luminosity functions of the SDSS DRl groups/clusters. Thin lines show the luminosity functions found using 
the clusters with at least two galaxies in the observational window, bold lines show luminosity functions for all groups/clusters, 
including groups with only one galaxy in the observational window. When calculating the total luminosities of groups, the 
Schechter function parameters of the set El were used. Right: the cumulative function of the total masses of DM-haloes of the 
model M200B. 



procedure (for details see E03a). The luminosity function was 
calculated separately for groups/clusters with at least two vis- 
ible galaxies, and for all groups/clusters including the systems 
with only one visible galaxy in the visibility window. In both 
cases the numbers of clusters have been corrected for selec- 
tion effects. We see that in the second case the number of 
groups/clusters per unit volume is larger by a factor of ~3 for 
the low luminosity section of the luminosity function. This re- 
sult shows that the groups with one bright main galaxy domi- 
nate among low-luminosity groups. We believe that this higher 
density represents the true number-density of low-luminosity 
groups better than the density found from groups with at least 
two galaxies in the observational window. It is well known that 
in our vicinity the majority of groups are similar to our Local 
Group, which consists of two subgroups with one bright main 
galaxy (our Galaxy and M31) and a number of considerably 
fainter companion galaxies. 



3.4. Luminosities of galaxies in different environments 

It is well known that in clusters brighter galaxies are concen- 
trated toward cluster centres. For this reason the smoothing 
scale must be relatively small in order to be sensitive to the po- 
sitions of galaxies within clusters. Taking this into account we 
shall use the density found with the 2 /; ' Mpc smoothing as an 
environmental parameter to describe the surrounding density of 
galaxies. The luminosity of galaxies as a function of the envi- 
ronmental density is shown in Fig. |5l here we used the SDSS 
EDR to calculate the luminosity density. This figure shows that 
the most luminous galaxies in high-density regions are about 
5-10 times more luminous than the most luminous galaxies in 
low-density environments. It is interesting that the contrast in 
luminosity between the high- and low-density regions is of the 
same order than the luminosity contrast of groups/clusters be- 
tween the high- and low-density regions. 
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Fig. 5. The luminosities of galaxies in various density environ- 
ments for the SDSS EDR Southern slice. The environmental 
density was calculated using Gaussian smoothing with the rms 
scale of 2 Mpc. 



4. Clusters in simulations 

4.1. DM lialoes and density fields 

The second goal of our study is comparison of observational 
data with numerical simulations. We have used for this purpose 
three simulations. Ml 00, M200A and M200B, with 128^ and 
256"' particles in 100 Mpc and 200 /i"' Mpc cubes, respec- 
tively. The cosmological parameters of the models are given in 
Table 12 In the analysis that follows we used the DM-haloes 
identified by the FoF algorithm with the search radius 0.23 (in 
units of the mean particle separation) for the model M200B. 

The density fields for the N-body simulations were found 
using all particles of simulations, i.e. we calculated the sim- 
ulated true total matter density fields. The high-resolution 
density was found using the conventional cloud-in-cell (CIC) 
scheme and additional Gaussian smoothing with the rms scale 
0.8 in grid cell units (0.6 li^^ Mpc). Gaussian smoothing was 
used in order to avoid the presence of numerous empty cells in 



6 



J. Einasto et al.: Environmental Effects in SDSS clusters 




Fig. 6. The dependence of the parameters of DM-haloes in the model M200B on the density of the environment, D (in units of 
the mean density). The density of the environment was determined from the low-resolution field using the Epanechnikov kernel 
with a radius of 8 Mpc. Left panel: total mass; right panel: full velocity. 



low-density regions. The low-resolution density field was cal- 
culated using an Epanechnikov kernel with the radius 10 in grid 
units, which corresponds to 8 /i ' Mpc. This field was applied 
to find simulated superclusters and environmental densities of 
DM-haloes. 

4.2. Properties of DM-haloes in different environments 

The dependence of the total mass of DM-haloes on the den- 
sity of the environment is shown in the left panel of Fig. |6l 
Here the dependence of the DM-halo mass on the density of 
the environment is very well expressed: the most massive clus- 
ters in a high-density environment are by a factor of a hundred 
more massive than the most massive clusters in a low-density 
environment. For a numerical simulation, we have also infor- 
mation on the velocities of DM-haloes, which is absent for the 
SDSS DRl groups/clusters. The right panel of Fig. |6l shows 
the dependence of the full cluster velocity on the density of the 
environment. The full velocity has a less pronounced density 
dependence, but here, too, DM-haloes in most dense environ- 
ments have a factor of ten larger velocities than DM-haloes in 
less dense environments. 

A similar result was obtained by E04b using a different def- 
inition of the density of the environment (the distance to the 5th 
nearest neighbour). In this paper, in addition to the correlations 
considered here, the virial radii and the intrinsic rms veloci- 
ties of DM-haloes in various environments were also studied. 
The virial radii of DM-haloes in a high-density environment 
are larger than those in a in low-density environment, but here 
the contrast between the high- and low-density regions is not so 
large. The rms velocities of DM-haloes have a very strong en- 
vironmental dependence, similar to the dependence observed 
for masses. This is natural, as the rms velocities of virialized 
haloes are determined by their masses. 

5. The evolution of various environmental regions 

To understand better the dependence of cluster properties on 
the environment we determined the distribution of particle den- 



sities in regions of various environmental density. For this pur- 
pose, for the model M200A we found for every particle in 
the simulation two density values, the local and the global 
density. The local density attributed to the particle was found 
as described above with Gaussian smoothing with the rms 
scale 0.8 in grid cell units (0.6 Mpc). The global density 
was found using the low-resolution density field as described 
above (smoothed with the Epanechnikov kernel of the radius 
8 /i ' Mpc). The density fields and particle densities were found 
for four epochs of the simulation, corresponding to the redshifts 
z = 5, 2, 1,0. The simulations started at the redshift z - 50, 
so at all epochs considered in our analysis the density field was 
well evolved. 

We divided the whole simulation volume into four re- 
gions by increasing global density. These regions correspond to 
voids, poor and rich filaments, and superclusters. The analysis 
of the distribution of particles for the present epoch z = shows 
that approximately 50% of all particles are presently located 
in rich supercluster regions with the global density D > 2.661. 
Rich filaments (actually poor superclusters) can be localised as 
systems lying in the range of global densities 1 < D < 2.661; 
at the present epoch about 25% of particles lie in this density 
range. Poor filaments can be found as systems in the global 
density range 0.484 < D < 1; about 15 % of particles lie 
presently in this density range. The rest of particles at global 
densities D < 0.484 constitute the void region; the fraction of 
particles in this density range is at the present epoch about 10%. 

If we are interested in the dynamics of different individ- 
ual regions we have to trace back the positions of individ- 
ual particles. This approach has been followed by Gottlober 
et al. (2003 1 in their study of the evolution of individual voids. 
Our task is simpler, as we are interested in the evolution of 
the simulation sample as a whole. It is clear that during the 
evolution DM-particles cluster locally to form DM-haloes, but 
most particles do not change their large-scale environment. The 
larger the scale, the smaller are the velocities at that scale. In 
other words, we may assume that the same fraction of particles 
presently located in the region of the highest global density was 
in earlier epochs also in the regions of the highest global den- 
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Sample 


Threshold 


Volume 


Density 






Do Oi D2 


Void Pfil Rfil Scl 


Void P fil R fil 


Scl 



z = 5 0.706 0.845 1.035 17.3 20.8 27.1 34.8 0.548 0.732 0.929 1.442 

z = 2 0.587 0.828 1.239 25.8 24.8 26.3 23.0 0.374 0.619 0.955 2.166 

z=l 0.527 0.857 1.552 32.8 26.9 24.4 15.9 0.289 0.566 1.035 3.151 

z = 0.484 1.000 2.661 46.7 26.7 19.1 7.5 0.207 0.549 1.350 6.641 



sity, and similarly the fractions of particles in other ranges of 
global density did not change. Under this assumption we found 
for each simulation epoch threshold global densities D, which 
divide the sample of all particles at a given epoch into regions 
of global density which occupy, from lowest values upward, 
10%, 15%, 25%, and 50% of all particles. As noted above, we 
call the respective regions the void regions, the poor and rich 
filament regions, and the supercluster regions. The respective 
global density thresholds for all epochs considered are given in 
TableOl 

During their dynamical evolution, superclusters shrink in 
volume and voids expand. To follow the change of the volume 
of various environmental regions we give in Table|3lthe volume 
occupied by various regions (in per cents of the total simulation 
volume). The distribution of volume fraction as a function of 
the local density d is shown in Fig. for all epochs of simu- 
lation. Different lines show the distribution for various global 
density D ranges: the void, the poor and rich filament, and the 
supercluster range. The fraction of particles in each region is 
known, so we can calculate the mean density of matter in each 
region (in units of the overall mean density); these mean densi- 
ties are also given in the table. We see that void regions occupy 
initially (at z = 5) only 17% of the volume, during evolution 
this fraction grows to 47%, and the mean density shrinks from 
0.5 to 0.2. Fig.0shows that in lowest density regions the den- 
sity is at present epoch less than 0.01. On the other hand, the 
50% of all matter that lies in high-density regions (superclus- 
ters) occupies initially 35% of space, this fraction decreases to 
8%, and respectively the mean density increases from 1.4 to 
6.6; in the highest density regions it is higher than 100. The 
poor and rich filament regions have intermediate behaviour, 
their volume fractions and mean densities changing much less. 

We determined in each of the regions the distribution of 
the particles by their local density d. Using the local density, 
we can assign particles to different populations. We emphasize 
that particles with local densities less than 1 cannot belong to 
clusters or groups, since galaxy formation starts only when the 
local density exceeds a certain critical threshold much higher 
than the mean density (see Press & Schechter ll974> . Following 
these ideas we classify the population of particles with local 
density below unity as primordial (non-clustered) particles, the 
population of particles with the local density values in the range 
1 <d< 10 as poor cluster (group) particles, and the population 
of particles with the local density values d > 10 as rich cluster 
particles. For an illustration of the distribution of particles of 



populations with various local densities see Fig. 1 of Einasto et 
al. (IT^ . 

The distribution of the number of particles according to 
their local density d is given in Fig.Ofor all simulation epochs 
considered. The Figure is similar to the previous one, only here 
we plot the distribution of mass instead of volume. The Tabled] 
gives the fractions of particles (in per cents of the total number 
of particles in simulation), which belong to populations of var- 
ious local density: the populations of the primordial, the poor 
cluster and the rich cluster particles; these populations are de- 
noted as 0, 1 and 2, respectively. The fractions are given, as 
explained above, separately for the void, the poor filament, the 
rich filament, and the supercluster environment. 

Let us discuss the data in more detail. Consider, first, the 
void regions. The void particle population distributions are 
drawn in Figs. and |S]by solid lines (see also the columns 
labeled in Table |^. We see that at all epochs the distribu- 
tion of local densities of void particles is rather symmetrical, 
both in volume and mass. In void regions most particles belong 
at all epochs to the primordial non-clustered population 0. At 
the epoch z - 5 about 1 % of the particles belong to the poor 
cluster population (the columns labeled 1 in Table|4}; this frac- 
tion increases with time and reaches 2.2% at the present epoch 
{z - 0). Initially there are no particles of rich cluster popula- 
tion among void particles (columns 2), and at the present epoch 
a very tiny fraction of particles has crossed the threshold of the 
rich cluster population. 

The distribution of particles in the poor filament regions is 
plotted in Figs.0and |8]by dashed lines. We see that the initial 
distribution of particles is similar to the distribution in void re- 
gions; however, the distributions are shifted toward higher local 
density values, and thus the fraction of the poor cluster popu- 
lation 1 is higher The poor cluster population 1 grows with 
time, so that at the present epoch about half of the poor filament 
particles belong to it. The rich cluster population 2 fraction in 
these regions grows from zero to 1.5%. Most of the volume in 
the poor filament region is occupied by local voids with local 
density values d <l (see Fig.0. 

The evolution of particles in the rich filament environment 
is shown in Figs.0and |Slby dot-dashed lines. In this region ini- 
tially about half of the particles belong to the primordial pop- 
ulation 0, and the other half to the poor cluster population 1 ; 
there are no particles of the rich cluster population 2. As time 
goes by, the fraction of primordial particles rapidly decreases 
and the fraction of rich cluster particles increases; the fraction 
of poor cluster particles changes little. 
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Fig. 7. The distribution of volume fraction as a function of local density d in regions of various environment in model M200A. The 
distribution is found for various regions of global density which correspond to voids, poor and rich filaments, and superclusters, 
containing 10%, 15%, 25%, and 50% of all particles in regions of increasing global density D, respectively. The upper panels are 
for the epochs z = 5, and z = 2, and the lower panels for the epochs z = 1 and z = (from left to right). 



Table 4. The fraction of particles of various local density in different environments 



Sample 





Voids 
1 


2 


Poor filaments 
1 2 


Rich filaments 
1 2 


Superclusters 
1 2 


z = 5 


8.4 


1.1 


0.0 


10.4 


4.8 


0.0 


11.6 


13.5 


0.0 


8.9 


40.1 1.1 


z = 2 


8.3 


1.4 


0.0 


8.8 


6.5 


0.1 


8.1 


16.1 


1.0 


4.8 


26.4 18.7 


z= 1 


7.9 


1.6 


0.0 


7.4 


7.4 


0.4 


5.8 


15.3 


4.2 


2.6 


16.8 30.7 


z = 


7.4 


2.2 


0.0 


5.0 


8.1 


1.5 


3.0 


12.9 


9.9 


0.8 


8.8 40.3 



For the evolution of the supercluster region see the dot- 
dash-dashed lines in Figs. and |8] In this environment the 
fraction of primordial particles rapidly decreases with time to 
almost zero, the fraction of poor cluster particles 1 decreases 
from about 40% to 9%, and the fraction of rich cluster popu- 
lation 2 increases from 1% to 40% (from z = 5 to z = 0). A 
large fraction of particles have local densities far in excess of 
our limiting density d - 10. However, a considerable fraction 
of space is still occupied by local voids, as seen from Fig.0 

These results can be summarized as follows: in void regions 
the mean density decreases continuously, as a result DM-haloes 
almost do not evolve dynamically, and most particles remain 
as primordial (non-clustered) ones. In supercluster regions the 



dynamical evolution is very rapid, the primordial population 
clusters rapidly, and later evolution consists of the transition 
of galaxies and groups to rich clusters. Here we have not fol- 
lowed the evolution of individual clusters, but it is clear that 
in this later stage merging of smaller DM-haloes to form rich 
DM-haloes plays an important role. In poor and rich filament 
regions the evolution is between these two extreme cases. 
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Fig. 8. The distribution of particles as a function of local density d in regions of various environment in the model M200A. 
The distribution has been found for various regions of global density which correspond to voids, poor and rich filaments, and 
superclusters, containing 10%, 15%, 25%, and 50% of all particles in regions of increasing global density D, respectively. The 
upper panels show the epochs z - 5, and z = 2, and the lower panels show the epochs z - I and z = (from left to right). 



6. Discussion 

6.1. Luminosity/mass functions in real and simulated 
clusters 

Let us compare first the luminosity and mass functions of real 
and simulated groups/clusters. In our preliminary study (E03a 
and E03b) we defined groups and clusters as enhancements 
of the 2-dimensional high-resolution luminosity density field. 
In the present work we used full 3 -dimensional data to define 
groups/clusters, both for the real data and simulations. 

New more accurate data show (see Fig. 0}, that there is 
no essential difference between the luminosity functions in the 
Northern and Southern strip of the survey. Our preliminary 
analysis based on the SDSS EDR and groups/clusters defined 
using the 2-dimensional luminosity density data suggested the 
presence of differences between the Northern and Southern 
strips. The new analysis does not support our previous conclu- 
sion. 

The volume density of groups/clusters according to the 
SDSS DRl data is 3 x 10"^ (/i ' Mpc)"^ for L > 10'' Lq 
groups/clusters. This estimate is in fairly good agreement with 
the estimates of the number density of groups based on the 
group mass function by Girardi & Giuricin (.2000 1, see also 
Heinamaki et al. (.2003 1. 



The luminosity of simulated DM-haloes is not well- 
defined. Thus we use for comparison the integrated mass func- 
tion of DM-haloes in simulations, presented in the right panel 
of Fig.|4] We see that the overall shape of the integrated mass 
function is rather similar to the luminosity function of the real 
groups/clusters. In the low-mass range the mass function of 
DM-haloes is steeper than the luminosity function of the real 
groups/clusters. The presently available data are insufficient to 
judge if this diff'erence is significant or not; a more detailed 
study of simulated samples is needed, with simulated galax- 
ies generated in DM-haloes. The mean volume density of DM- 
haloes is very close to the mean volume density of groups in 
the SDSS DRl, thus we can say that the population of DM- 
haloes represents the population of real groups/clusters rather 
well (the overall bias is small). 

6.2. Environmental effects in real and simulated 
clusters 

Let us compare now the environmental dependence in real and 
simulated cluster samples. The environmental dependence has 
been investigated using two completely independent param- 
eters to characterize the large-scale environment. E03c com- 
pared groups and clusters in high density environments, de- 
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fined as a neighbourhood of rich clusters, and in low density 
environments, far from rich clusters. E03d compared the prop- 
erties of groups and clusters that belong to superclusters and 
the properties of groups/clusters that do not belong to super- 
clusters. E04b used the distance to the 5th nearest neighbour as 
a parameter of the large-scale environment to describe the en- 
vironment of DM-haloes in simulations. E03a and E03b used 
for this purpose the low-resolution luminosity density field, as 
we do in the present study. 

These parameters of the large-scale environment used in 
various studies are independent of each other and characterize 
the environment from different points of view. The luminosity 
density field approach takes into account the luminosity or the 
mass of neighbouring objects, including the dark matter parti- 
cles in simulations. The nearest neighbour approach, as well as 
the proximity to rich clusters, depends only on the number and 
position of neighbours. Their luminosity or mass is ignored. 

All relations considered so far between the physical pa- 
rameters of groups/clusters and the environmental parameter 
show the presence of well-defined correlations, both in real as 
well as in simulated samples: in high-density regions (super- 
clusters) groups/clusters are brighter and more massive than 
in low-density regions (void regions). Also they have slightly 
larger radii and greater bulk velocities (see E04b). 

6.3. Comparison to previous work 

Historically, the dependence of galaxy properties on their large- 
scale environment has been investigated long ago, starting from 
the pioneering studies by Davis & Geller ( I1976> and Dressier 
J198 0i. In these early studies a striking contrast between the 
morphological types of galaxies in the cluster and field envi- 
ronments was found: in clusters elliptical galaxies dominate 
and in the field dominate spiral and irregular ones. More re- 
cently Einasto (199 la 1991b) and Mo et al. ( 1992i found the 
dependence of galaxy luminosity and morphological type on 
the large-scale environment up to the scale of 10 /z ' Mpc. This 
result was confirmed by Lindner et al. (1995 1996) by the 
study of voids defined by galaxies of different absolute magni- 
tude: bright galaxies define much larger voids than intrinsically 
faint galaxies; these faint galaxies form poor filaments inside 
large voids surrounded by bright galaxies. 

Recently Balogh et al. (2004), and Blanton et al. (2004a 
12004b, see also references in these papers) investigated the 
dependence of physical properties of galaxies in different lo- 
cal and global environment. Among physical properties they 
considered colours, luminosities. Ha emission, and the Sersic 
<1968> luminosity radial profile index. The local environment 
was defined as the spatial density on the 0.5-1 Mpc scale, 
the global density on the 5-10 /; ' Mpc scale. Blanton et al. 
J2004bt come to the conclusion that the blue galaxy fraction, 
and the recent star formation history in general, depend mostly 
on the local environment. On the other hand, Croton et al. 
J2004> find that the luminosity function of galaxies depends 
strongly on the global environment (see below). 

Peebles (2001) compared void and cluster galaxies and 
confirmed the presence of a striking difference of the proper- 



ties of these galaxies. He argued that this difference may be a 
challenge to the ACDM model of structure formation. Mo et al. 
("20041 came to the conclusion that the differences in luminos- 
ity observed in regions of different environmental density can 
be attributed to the differences in the mass of DM-haloes from 
which galaxies form. As our present study shows, the masses 
themselves depend on the environment. 

The luminosity dependence of galaxy clustering was inves- 
tigated by Norberg et al. (12001 1 i2002i) . using the 2dF Galaxy 
Redshift Survey. They found that the clustering amplitude in- 
creases with absolute magnitude, confirming earlier results by 
Einasto, Klypin & Saar (IT986b . Bromley et al. ( [T998b . and 
Beisbart & Kerscher (2000). 

The luminosity function of galaxies in different environ- 
ment was recently investigated by Hiitsi et al. ( I2002I I. De 
Propris et al. 12003 1 and Mo et al. (©041. Hutsi et al. esti- 
mated the luminosity function in three regions of global den- 
sity, defined by the \Q h^^ Mpc smoothing, and found that the 
characteristic luminosity of galaxies increases by a factor of 
1.5 from the low-density to the high-density regions. Mo et al. 
argued that this difference in the characteristic luminosity is in 
agreement with their model prediction, which assumes that the 
Schechter approximation is valid in virtually all environments. 

A very detailed study of luminosity functions of galaxies in 
regions of different density of the large-scale environment was 
made by Croton et al. ( 2004 1, using the full dataset of the 2dF 
Galaxy Redshift Survey. Using densities smoothed on a scale 
of 8 /i ' Mpc as in the present paper, he divided the volume 
under study into 7 regions of various density of the environ- 
ment, from extreme void to cluster populations. In all regions 
the luminosity function was calculated and the parameters of 
the Schechter function were found. The bright end of the func- 
tion depends primarily on the characteristic absolute magnitude 
M* of the Schechter function. For the extreme void population 
this parameter is -18.3, and for the cluster population -20.1. In 
other words, the brightest galaxies in voids are approximately 5 
times fainter than in clusters. This result is in very good agree- 
ment with our data on the distribution of galaxy luminosities of 
the SDSS in various environments (see Fig.|5j). 

The dependence of the total luminosities and masses of 
galaxy systems on the density of the environment has been in- 
vestigated only recently by E03a, E03b, E03c, E03d, E04a and 
E04b (for a discussion about the properties of haloes in dif- 
ferent environments see references in this paper). These stud- 
ies show a tendency similar to that of galaxy luminosities - in 
high-density regions massive and luminous clusters dominate, 
whereas in low-density regions all galaxy systems are poor and 
faint. The most luminous clusters in a high-density environ- 
ment are a factor of 5-10 more luminous than the most lumi- 
nous clusters in a low-density environment. It is striking that 
this factor is in the same range as for the most luminous galax- 
ies in high- and low-density environments. 

The properties of groups of galaxies in the vicinity of rich 
clusters were compared with the properties of ordinary groups 
by Ragone et al. ( I2004t . They used a sample of groups iden- 
tified in the 2dF Galaxy Redshift Survey, and compared the 
properties of groups in the vicinity of rich clusters and the rest 
of the sample. The observational results were compared with 
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simulated clusters using the Virgo Consortium Simulation. In 
both the real and simulated groups there exist similar relations: 
the larger the host mass, the higher is the luminosity or the mass 
of the DM-halo. 

Comparison of the richness of DM-haloes in different en- 
vironments was made by Gottlober et al. (2003 1, using high- 
resolution numerical simulations. Their results show that DM- 
haloes in voids have much lower masses than in high-density 
environments. This result was recently confirmed by Colberg 
et al. (2004 ). Their Fig. 14 demonstrates that the most massive 
DM-haloes in voids are about 100 times less massive than the 
most massive DM-haloes in general; also the growth of masses 
of DM-haloes in voids is less rapid than in general. 

6.4. Interpretation of the environmental dependence of 
galaxy and cluster properties 

One may ask, why are void galaxies and clusters so different 
from galaxies and clusters in dense regions? 

The evolution of DM-haloes in numerical sim- 
ulations has been investigated by a number of au- 
thors. Several of these simulations have been vi- 
sualized, as an example those by Andrey Kravtsov 
(http : //cfcp . uchicago . edu/lss/f i laments .html i. 
In this simulation the formation of rich clusters and super- 
clusters consisting of a system of intertwined filaments can be 
clearly seen. 

We are interested in the difference between the structure 
of galaxy systems in high- and low-density environments. This 
problem has been studied in detail by Gottlober et al. ( 2003, 
for earlier work see references in this paper). They chose sev- 
eral large under-dense regions of a diameter of ~20 /z ' Mpc 
(voids defined by bright simulated galaxies) and re-simulated 
the evolution of these voids with a very high mass resolution 
4.0 X 10^ /z^'Mq. Their Fig. 2 shows the distribution of dark 
matter in one of these voids at the epochs z - 2 and z = 0. At 
both epochs a well-developed system of filaments with com- 
pact knots (DM-haloes) can be seen. However, the masses of 
these haloes are at both epochs of the order of 10'^ 'Mq, only 
the most massive ones have masses of few times of 10'" /i^'Mq. 
In other words, after the early growth the knots stop growing. 
The volume density of DM-haloes in these voids is a factor of 
10 lower that in the simulation as a whole. 

In contrast, the growth of DM-haloes in high-density re- 
gions is much more rapid and continues over the whole period 
under study. As shown, among others, by Frisch et al. (|1995 1, 
the first objects to form in simulations are rich clusters in su- 
perclusters. Our calculations presented in the last Section fully 
confirm these earlier conclusions. 



7. Conclusions 

The main results of our study of clusters and superclusters in 
the SDSS DRl and the comparison with the results of numeri- 
cal simulations can be summarised as follows: 



- We have found groups and clusters in the SDSS DRl data 
using three-dimensional information on the distribution of 
galaxies. 

- Using Gaussian smoothing with the rms scales of 0.8 and 
10 /z ' Mpc we have derived high- and low -resolution lumi- 
nosity density fields for the SDSS DRl data in two equa- 
torial strips; the low-resolution density was used as an en- 
vironmental parameter to describe the large-scale environ- 
ment of groups and clusters. 

- New three-dimensional data confirm the environmental de- 
pendence found earlier using the two-dimensional data: in 
high-density regions (superclusters) groups and clusters are 
richer and bigger, and galaxies themselves are brighter. 

- Numerical simulations show a similar tendency: in a high- 
density environment DM-haloes are richer, and have larger 
bulk motions than DM-haloes in a low-density environ- 
ment. 

- Our explanation of the density-luminosity relationship is 
by the combined influence of density perturbations of all 
scales. Superclusters are the regions where the density per- 
turbations of large and medium wavelength combine to 
generate high-density peaks: here the overall density is 
high and the dynamical evolution of clusters and galax- 
ies is rapid and continues until the present. Voids are re- 
gions where large-scale density perturbations have negative 
amplitudes; here, due to medium and small-scale perturba- 
tions also a filamentary web forms; however, due to the low 
mean density the dynamical evolution is slow and stops at 
an early epoch. 
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